Scaling of dynamics with the range of interaction in short-range attractive colloids 
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We numerically study the dependence of the dynamics on the range of interaction A for the short- 
range square well potential. We find that, for small A, dynamics scale exactly in the same way as 
thermodynamics, both for Newtonian and Brownian microscopic dynamics. For interaction ranges 
from a few percent down to the Baxter limit, the relative location of the attractive glass line and 
the liquid-gas line does not depend on A. This proves that in this class of potentials, disordered 
arrested states (gels) can be generated only as a result of a kinetically arrested phase separation. 

PACS numbers: 61.20.Ja, 82.70.Dd, 82.70.Gg, 64.70.Pf 
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Colloidal dispersions form gels, disordered arrested 
states of matter at low packing fraction (j>, if the colloid- 
colloid hard sphere repulsion is complemented by a short- 
range attraction iL, ^ . The nature of the gel tran- 
sition in short-range attractive colloidal systems has re- 
ceived significant attention in recent years (for a recent 
review see for example Ref. j5j). Several routes to the 
gel state have been proposed and critically examined. In 
particular, it has been speculated that the gel-line consti- 
tutes the extension to low cj) of the attractive-glass line, 
an idea which would provide an unifying interpretation 
of the gel and glass arrest state of matter |^ . An alter- 
native interpretation suggests that colloidal gel results 
from an interrupted liquid-gas phase separation, inter- 
rupted by the glass transition which takes place in the 
dense regions created during the spinodal decomposition 
kinetics l3illiP|- The two scenarios, which differ only by 
the relative location of the glass line(s) with respect to 
the phase separation line, are sketched in FigQ In case 
(i), the attractive glass line pre-empts the meta-stable 
liquid-gas separation and the gel line can be approached 
from equilibrium conditions (Fig^). In case (ii), the 
glass line meets the phase separation line on the high cj) 
branch, and the morphology of the low (j) arrested state 
is dictated by the phase separation kinetic (Fig^). 
The thermodynamic phase diagram of simple models for 
short-range attraction has been evaluated theoretic ally 
and successfully compared with experimental data |l3l 
^4, 15, 16] . When the attraction range is a few percent 
of the particle size, the equilibrium phase diagram is com- 
posed only by a fluid-phase and a crystalline phase. The 
liquid-gas coexistence locus is hidden within the region of 
fluid-crystal coexistence. For small range of attractions, 
the liquid-gas coexistence curve for different models can 
be scaled onto each other by comparing different systems 
at the same value of the second virial coefficient , pro- 
viding an effective characterization of the dependence of 
the liquid-gas coexistence line on the range of attraction. 

The dependence of dynamic properties on the range of 
attraction has been studied at large (j) within the mode- 
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FIG. 1: Sketch of two possible relative location in the T-(j) 
plane of the liquid-gas coexistence line and of the glass lines. 
In a) the liquid-gas coexistence is hidden below the liquid- 
glass transition line (case (i) in the text). In b) the glass line 
intersect the binodal below the critical temperature at a (/> 
value higher then the critical one (case (ii) in the text). The 
As point is the MCT singularity 111 III that might be 
encountered for the range values we discuss in this work. 



coupling theory For small ranges, two distinct glass 
lines appear, indicating the possibility of forming two 
distinct glass states, commorily named repulsive and at- 
tractive glasses m El El IP . According to 
MCT lOjlL 12»,2M^ for small attraction range, the at- 
tractive glass line extends to low 0, supporting type (i) 
scenario. Unfortunately, MCT overestimates the loca- 
tion of the glass lines and a non negligible mapping must 
be applied on theoretical curves before com par ing the- 
ory with experimental or simulation data |23, 123 ■ In the 
case of a square well potential with an attractive range 
of 3% of the hard-sphere diameter such a mapping has 
been evaluated and the mapped attractive glass line has 
been found to end on the right side of the spinodal , 
in agreement with type (ii) scenario. 
If the location of the attractive glass line and of the 
liquid-gas line depends in different ways on the range of 
the attraction a transition from case (ii) (which is known 
to be the correct case for interaction ranges larger than 
3% ,9| ) to case (i) could take place at a very small crit- 
ical value of the attraction range. This Letter addresses 
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this question, by examining the dependence of the dy- 
namics on the attraction range, both for Brownian (BD) 
and Newtonian (ND) dynamics. We show that, for in- 
teraction ranges from a few percent down to the Baxter 
hmit (2Jj and for packing fraction smaller than 0.40, (the 
packing fraction range where gels are observed in exper- 
iments), dynamics and thermodynamics loci scale with 
the range of interaction in the same way, ruling out case 
(i) as route to the gel formation in short-range attractive 
colloids. 

We investigate a system that has been extensively 
studied earlier, a binary square well (SW) mixture |2^ 
|2^ . The binary system is a 50%-50% mixture of = 
2000 particles. The two species (labeled A and B) are 
characterized by a diameter ratio uaIob = 1.2. Masses 
are chosen to be equal and unitary. The attraction is 
modeled by a SW interaction defined according to: 

{oo r < (T„,/3 
-uq ao,,f} <r <ao,,s3 + l^cS (1) 
r > (Ja.jS + ^a,l3 

where Ua^p — [oa +o'/3)/2, a, (3 — A, B and ^a,p is the 
range of the attraction. We fix ctq^^ and vary the relative 
well- width e = ^ "^^^ — -. We report data for extremely 

small relative well width — from 10~^ to 5 • 10~^ — 
covering an interval which starts from a physically re- 
alizable limit and ends close to the theoretical Baxter 
limit. We choose /cs = 1 and set the depth of the po- 
tential uq = 1. Hence T = 1 corresponds to a thermal 
energy fc^T equal to the attractive well depth. The di- 
ameter of the small specie is chosen as unity of length, i.e. 
ub = 1. Density is expressed in term of packing fraction 
4> = {pA<j\ + Pb(J%) ■ 7r/6, where pa = N^/L^, L being 
the box size and Na the number of particles of species 
a. Time is measured in units of ctb • (m/uo)^^'^ . ND 
has been coded via a standard event driven algorithm, 
commonly used for particles interacting with step-wise 
potentials BD has been implemented via the posi- 
tion Langevin equation: 

Mt)-^m+rdt), (2) 

coding the algorithm developed by Strating . In Eq[21 
ri{t) is the position of particle i, ii{t) is the total force 
acting on the particle, Dq is the short-time (bare) diffu- 
sion coefficient, ri{t) a random thermal noise satisfying 

< ri(t)ri{0) >— kBT5{t). In Strating's algorithm, a ran- 
dom velocity (extracted from a Gaussian distribution of 
variance ksT /m) is assigned to each particle and the 
system is propagated for a finite time-step , accord- 
ing to event-driven dynamics. We chose such that 
short time motion is diffusive over distances smaller than 
the well width. For the smallest e, reliable estimates of 
dynamical properties require more than lO"'^" collisions 
(about two weeks on a 3GHz processor). 

For the very small values of relative well width e con- 
sidered here, thermodynamic properties at different well 



width can be scaled by using as scaling variable the value 
of the second virial coefhcient B2 |3 • For the SW binary 

B^^ 4-B^ ^ +2B^^ 

mixture B2 = — \ — where 

B^'^ = l^<p [1 - (e^"° - - e)-' - 1]] • (3) 

For the equivalent 50-50 hard-sphere binary mixture, 
B2 is 

j^HS ^ I Ha+^BB+^^bI 

An adimensional second virial coefhcient can be defined 
as B2 = B2/B2^ . This quantity helps in comparing 
between different models and different samples [23 • At 
small e, B2 becomes essentially function of the variable 
gg/3uo Ae'^"" . In the same limit, state points at the 
same B2 and cj) are characterized, to a very good approx- 
imation, by same thermodynamic properties, i.e. same 
bonding pattern, same energy, same structure. In the 
limit e — > the system behaves similarly to the Bax- 
ter model [3| at the same B2 state point. The Baxter 
potential VB{r) is best defined via 

12r 

where t is the adhesiveness parameter, which plays the 
role of effective temperature, 9 and 5 are respectively 
the Heaviside and Dirac functions. This model has been 
extensively used in the interpretation of experimental 
data [s^l despite its known pathologies 1^41 ■ For the 
Baxter potential, B2 — 1 — 1 /4r and the location of the 
liquid-gas critical point, recently determined with great 
accuracy, is (f)c = 0.266 and B^ = -1.2 

Fig. 121 shows the spinodal line for the SW model with 
e = 0.01, estimated by bracketing it with the lowest T 
stable point and the first phase separating state point 
along each isochore. It also show the data from Miller 
and Frenkel for the Baxter potential. The agreement 
between the two set of data, notwithstanding the differ- 
ent system (binary mixture vs. one component, SW vs. 
Baxter) confirms that the Baxter limit is already reached 
when e = 0.01. Fig. |21 also shows the ND isodiffusivity 
lines [s^l, defined as the locus where the normalized dif- 
fusion coefficient D/Dq is constant. For ND, the nor- 
malization factor Dq = ^y^Sk^To^Jm = VthO'B accounts 
for differences in the microscopic time due to different 
thermal velocity Vt/^. The large values of D/Dq, even 
for « (/)c, confirms that, as in the previously studied 
e = 0.03 case 9], no arrested states can be approached 
in equilibrium for (j) < (pc- 

We next address the question of the dependence of the 
dynamics on e. We focus on two specific values of 0, 
respectively on on the left {(f) — 0.2) and on the right 
(0 = 0.4) of the critical point. For each </>, we select sev- 
eral pairs of e — T values such that B2 = —0.405. The 
average energy per particle are respectively of —2 and —4. 
Within our numerical precision, simulations for different 
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FIG. 2: Phase diagram for e = 0.01. The continuous line 
reproduces the coexistence curve calculated by Miller and 
Frenkel [s^ . The error bars represent our results, i.e. the 
intervals which bracket the spinodal (see text). Isodiffusiv- 
ity lines for three typical liquid values of D/Dq are plotted. 
Note that slow dynamics pre-glassy features (like two-step re- 
laxation decay) appear only when D/Dq < 5 ■ 10~*. 



e — T values converge to the same average potential en- 
ergy and same structure, supporting the hypothesis that 
for these small e values equality in B2 implies equal ther- 
modynamic properties. 

We focus on two dynamic quantities, the tagged- 
particle mean square displacement < r^(i) > and the 
bond autocorrelation function C{t), defined by: 
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FIG. 3: Mean square displacement for different e values at 
the same B| = -0.495 for ND (a) and BD (b) for cj> = 0.20 
and (p — 0.40. Do = vthO'B for ND and the bare monomer 
diffusion constant for BD. 
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where the N x N matrix Cij {t) defines the bonds at time 
t according to: 

. ^ _ J 1 if i and j particles are bonded 
c»jW-|o else 

Two particles are considered bonded if their relative dis- 
tance is in the attractive well. 

Fig.Elshows < r^{t) > for both ND and BD. Data are 
reported as a function of IDq to account for the trivial 
differences in vth for ND and in the bare self-diffusion 
coefficient Dq for BD. For both </> and both microscopic 
dynamics, < r^(t) > is independent on the range of the 
attractive potential, when the comparison is done at con- 
stant i?2- In other words, the only difference in the dy- 
namics is accounted by the trivial microscopic Dq scaling 
factor. This implies that the isodiffusivity curves calcu- 
lated for the e = 0.01 case, when reported in a i?^ — 
plane, describe the entire class of SW potentials with 
range shorter than e = 0.01. 

Fig.^shows C(i) as a function of IDq for different pairs 
(e,T) at fixed i?2- Since B2 is constant by construction, 
the average number of bonds in the system is the same for 
all investigated (e,T) pairs. In agreement with the data 
shown in Fig. O all C{t) collapse onto the same curve 



both for ND and BD. This suggests that, in IDq units, 
the probability of breaking a bond does not change along 
constant paths. It is worth stressing that, while the 
collapse is observed for both type of microscopic dynam- 
ics, the shape of the ND and BD correlation functions 
differs. In ND, C(t) is to a good approximation expo- 
nential while in BD it is stretched, with a stretching ex- 
ponent w 0.5. The same considerations hold for (j> = 0.40 
(Fig. ^p) . The fact that the decay of C [t] is still strongly 
affected by the microscopic dynamics, implies that MCT 
can not be applied at these (jj. 

In general, assuming that bond breaking is an acti- 
vated process, the bond breaking probability can be ex- 
pressed as a product of a frequency of bond breaking 
attempts uj times e~''"°, which express the probability of 
overcoming the barrier. In the case of ND, is pro- 
portional to the time requested to travel a distance of the 
order of A, and A/vth- Hence, the bond lifetime, 

apart from the thermal contribution vth absorbed in Dq, 
is controlled by the product Ae^"° , the same quantity 
controlling the value of B2 at small e and T. 
Results presented in Figs. EEl suggest that, for small e, 
the value of B2 characterizes not only thermodynamics, 
but also dynamics. In other words, at a given value of B2 



4 



U 



1 

0.8 
0.6 
0.4 
0.2 

1 

0.8 
0.6 
0.4 
0.2 





BD "\ \ND 

















...1 1 .1^ 



b) 

£=5x10 ' 
^ £=5x10"'^ 
£=5x10"" 
\-- £=5x10"^ 
xlO"^ 



10 



IQ" 10' 



\\ 



10 



10 



10 
tD. 



10 



10 



FIG. 4: Bond correlation function C as a function of tDo for 
different relative well width e and T but all at B2 = —0.405. 
Data for both ND and BD are reported, (a) (f) = 0.20; (b) 
(j) = 0.40. The insets show, for the case of ND, C{t) vs t. 



it corresponds a family of systems with different T and 
e, including the limiting case of Baxter, that posses the 
same static and dynamic properties. According to the 
present results, the apparent long bond lifetime charac- 
teristic of the Baxter model 33] is only induced by the 
extremely small thermal velocity associated to the van- 
ishing of T implicit in the limit e — ^ at fixed B2 ■ 
The similar scaling of dynamics and thermodynamics has 
important consequences for understanding gel formation 
in short-range attractive colloidal dispersions. The isodif- 
fusivity lines reported in FigI21 describe not only the case 
e = 0.01 for which they have been calculated but also 
the dynamics of all shorter ranged potentials, down to 
the Baxter hmit, at least up to the tested — 0.4 value. 
This has a profound consequence for the two scenarios 
discussed in Fig^ since it proves that in short-ranged 
potentials the glass line always meet the liquid-gas line 
on its right side. In this class of potentials, disordered 
arrested states at low (j) can only be created under out 
of equilibrium conditions, requiring a preliminary sep- 
aration into colloid rich (liquid) and colloid poor (gas) 
phases followed by an attractive-glass dynamic arrest in 
the denser regions. 
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